Predicting d° magnetism 
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Predicting magnetism originating from 2p orbit als is a delicate problem, which depends on the 
subtle interplay between covalency and Hund's coupling. Calculations based on density functional 
theory and the local spin density approximation fail in two remarkably different ways. On the one 
hand the excessive derealization of spin-polarized holes leads to half-metallic ground states and the 
expectation of room temperature ferromagnetism. On the other hand, in some cases a magnetic 
ground state may not be predicted at all. We demonstrate that a simple self-interaction correction 
scheme modifies both these situations via an enhanced localization of the holes responsible for the 
magnetism and possibly Jahn- Teller distortion. In both cases the ground state becomes insulating 
and the magnetic coupling between the impurities weak. 
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Conventional magnetism is associated with a narrow 
region of the periodic table, namely that of the 3d and 
5/ series. In these atomic shells strong Hund's coupling 
results in a high-spin configuration and therefore in the 
formation of localized magnetic moments. Then one also 
needs a mechanism for coupling those local moments. 
This is not universal and generally depends on the details 
of the compound investigated. Importantly the forma- 
tion of magnetic moments does not imply a macroscopic 
magnetic state. Furthermore, the orbital composition of 
the moment is crucial for establishing the microscopic 
nature of the coupling mechanism. 

Following the criterion of having large Hund's cou- 
pling one should expect to find high-spin configurations 
and possibly magnetism also in materials with partially 
filled 2p orbitals. This situation takes the name of d° 
magnetism 1 . Oxygen is the prototypical case: O2 has a 
triplet ground state and solid O orders in a complex an- 
tiferromagnetic structure at low temperature 2 . Another 
remarkable example is that of RD4O6, which was pre- 
dicted to be an anionogenic half- metal 3 , and measured 
spin-glass^. Away from molecular solids however the sit- 
uation is less clear. Two properties of standard solids 
conspire against p-type magnetism. First, the p shells are 
usually fully filled and therefore cannot sustain a mag- 
netic moment. Secondly the 2p hopping integrals are 
usually large and so are the typical bandwidths, mean- 
ing that the chances of moment formation are reduced 
by the large kinetic energy. 

Despite these complications the experimental claims 
for d° ferromagnetism are numerous and include 
graphitic carbon^ and both defective 6 - and doped oxides^. 
An explanation of these findings however remains con- 
troversial to date. The formation of the magnetic mo- 
ment is usually attributed to holes localized at the defect 
site, either this being the molecular orbital associated 
to cation vacancies^ or p-dopants at the site o 10 ! 11 . 
Then the magnetic coupling is justified with the degree 
of derealization of the same 2p shell responsible for the 
moment, in a scheme similar to the Zener model for stan- 
dard magnetic semiconductors^ 2 -. Critically most of the 
predictions are based on density functional theory (DFT) 
using local approximations of the exchange correlation 



potential (LSDA or GGA). These notoriously underesti- 
mate Coulomb repulsion and tend to over-delocalize the 
charge density. It is then no surprise that most of the 
calculations return a metallic (half-metallic) ground state 
and usually extremely large magnetic interaction. 

In this letter we demonstrate that, when strong cor- 
relation is introduced at the level of self-interaction cor- 
rection (SIC) the picture changes drastically and even 
the formation of a local moment becomes a delicate is- 
sue. Consequently the mechanism for magnetic interac- 
tion can be completely different from that described by 
the LSDA. We prove this idea by presenting three pro- 
totypical cases. First we look at F centers in Si02 for 
which LSDA predicts a delocalized hole sustaining mag- 
netic interaction instead of a paramagnetic self-trapped 
cente r 13 ! 14 . Then we look at B, C and N substituting for 
O in MgO, for which LSDA correctly predicts a magnetic 
ground state but fails in capturing the orbital order and 
the insulating state. Finally we look at the case of Ga 
vacancy and dopant Zn in GaN. The former case is simi- 
lar to that of F centers in Si02 while in the latter, LSDA 
now fails in predicting a high-spin state. 

Our calculations are performed using the standard 
DFT code siesta with LSDA functional^ and a de- 
velopment version implementing the atomic SIC scheme 
(ASIC) 16 . We treat core electrons with norm-conserving 
Troullier-Martin pseudopotentials, while the valence 
charge density and all the operators are expanded over 
a numerical orbital basis set, including multiple-^ and 
polarized functions^. The real space grid has an equiv- 
alent cutoff larger than 500 Ry. Calculations are per- 
formed with supercells of various sizes including /c-point 
sampling over at least 10 points in the Brillouin zone. 
Relaxations are performed with standard conjugate gra- 
dients until the forces are smaller then 0.04 eV/A. 

We begin our analysis by considering the famous pro- 
totypical case of the F center in Si02, i.e. Al substi- 
tuting Si (Alsi). In table [J we summarize the results of 
our LSDA and ASIC calculations. These have been ob- 
tained with a 36 atom unit cell containing a single Al 
impurity. The LSDA local geometry of Alsi has the four 
Al-0 bonds of equal length, 1.73 A, and the 5=1/2 hole 
spreads uniformly over the four O atoms. However relax- 



System cZlsda 



Mlsda 



^ASIC 



Masic 



Si0 2 : Al Si 1.73 (4) 1 

MgO: B 2.19 (6) 3 

MgO: Co 2.18 (2), 2.15 (4) 2 

MgO: No 2.10 (6) 1 
GaN: V Ga 2.10 (4) 
GaN: Znca 2.0 (4) 



1.67 (3), 1.94 (1) 1 

2.20 (6) 3 

2.18 (4), 2.15 (2) 2 

2.15 (4), 2.10 (2) 1 

3 2.14 (2), 2.28, 2.31 3 

1.96 (3), 2.4 (1) 1 



TABLE I: Summary table of the LSDA and ASIC calculated 
bond lengths, d (in A), and magnetic moments, \i (in /xb), for 
the various defects investigated in SiCb, MgO and GaN. In 
bracket the number of bonds of a given length. 



ation with ASIC produces a distortion. Three of the Al-0 
bonds relax to 1.67 A, while the forth gets considerably 
longer (1.94 A). Such a distortion is associated with the 
localization of the Al-induced hole, which occupies the 
orbital along the elongated bond. This picture is consis- 
tent with the experimentally found hole self-trapping and 
demonstrates the importance of SIC in the self-trapping 
problem^. 

If one now calculates the magnetic interaction between 
two Alsi, a remarkable result is found. By using a 72 
atom supercell containing two Al impurities placed at 
second nearest neighbours, we have calculated the en- 
ergy difference, AE, between the ferromagnetic and an- 
tiferromagnetic configurations of the cell. This is about 
120 meV for LSDA meaning that for Al concentrations 
around 8 % LSDA predicts Si02:Al to be a ferromag- 
net with an estimated mean-field critical temperature of 
240 K. No magnetism has ever been reported for Si02:Al 
at any concentrations and we conclude that this is sim- 
ply an artifact of the excessive LSDA hole derealization. 
When the same calculation is repeated with ASIC we find 
that the total energy difference drops to 1 meV, i.e. it is 
consistent with no magnetic coupling. These results set 
the theoretical benchmark of our method. 

Next we analyze the case of O substitutional defects 
A (A = B, C, N), in rock-salt MgO. These introduce 
respectively 3, 2 and 1 hole, so that the filling of the Ao 
2p shell is respectively 1/2, 2/3 and 5/6. In figure [1] we 
report the density of states for a 96 atom supercell doped 
with a single Ao and projected over the Ao 2p orbit als. 
For all the impurities both the LSDA and ASIC return a 
spin-polarized ground-state with a magnetic moment as 
large as the number of holes. The details of the charge 
distribution around Ao are however profoundly different 
in LSDA and ASIC. 

Let us look at the single acceptor No first. A magnetic 
ground state means that the hole is fully spin-polarized. 
The problem is then that of distributing a single hole 
among the three 2p-orbitals forming the chemical bonds 
with the Mg cations. As for Si02 also in this case LSDA 
predicts a spin-polarized ground state with the Fermi 
level cutting across a 2/3 filled minority N-2p impurity 
band. LSDA atomic relaxation is isotropic with a N-Mg 
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FIG. 1: (Color on line) Total density of states (black line) 
and density of states projected over the 2p shells of the anion 
acceptor (red line) for MgO doped respectively with B (a), 
C (b) and N (c). The left-hand (right-hand) side plots cor- 
respond to LSDA (ASIC). The majority (minority) DOS is 
plotted in the upper (lower) half of each panel. 



bond- length slightly larger (2. 10 A) than that between O 
and Mg of 2. 08 A. This half- metallic ground state, to- 
gether with the considerable energy overlap between the 
O 2p and N 2p orbitals (see FigJT]), is suggestive of a N 
ferromagnetic order via an impurity band Zener mecha- 
nism. In fact the calculated AE for a 96 atom cell and 
second nearest neighbours is 120 meV and ferromagnetic. 

Note however that this partially filled degenerate p 
configuration is sensitive to Jahn- Teller distortion. This 
is not captured by LSDA, due to the erroneous over- 
delocalization of the hole, but it can be described by an 
ASIC calculation. ASIC in fact leads to the expansion of 
the two of six N-Mg bonds (2.15 A) with the consequent 
N-2p levels splitting into a doubly degenerate occupied 
level just above the MgO valence band and an empty 
singlet. These are separated by a crystal field energy of 
approximately 3 eV. In this distorted configuration the 
magnetic moment is entirely localized over the longer of 
the p bonds as it can be clearly seen from Fig. [2] where 
we present the ASIC magnetization isosurface. 

As it is, MgO:N is an insulator and the only mecha- 
nism left for magnetic coupling between No is superex- 
change. This is expected to be extremely weak, since 
the first available Mg filled shells to mediate the virtual 
transition necessary for superexchange are the 2p, which 
are extremely deep, leading practically to no coupling. 
In fact the calculated coupling at second nearest neigh- 
bour is only 1 meV. Note also that, at variance to the 
standard Jahn- Teller distortion for d orbitals, the p-type 
distortion investigated here is frustrated, and each of the 
three degenerate 2p orbitals can transform into the sin- 
glet. This introduces an additional complication to the 






FIG. 2: (Color on line) Magnetization isosurfaces calculated 
with ASIC for No in MgO. Note that the 5=1/2 hole is com- 
pletely localized along the p y orbital and it is not evenly dis- 
tributed along all the bonds. Color code: blue=Mg, red=0, 
N is not displayed for clarity. 



superexchange mechanism, since only orbitals sharing the 
same angular momentum about the bonding axis can par- 
ticipate to the coupling. Finally we note that a similar 
Jahn- Teller distortion was predicted by Elfimov et al^ 
for N-doped SrO. 

Let us now move to MgO:C. This time the minority 
C-2p shell is singly occupied leading to a magnetic mo- 
ment of 2 /i B for both LSD A and ASIC. For Co also 
LSD A predicts Jahn- Teller distortion with four 2.15 A 
and two 2.18 A long Mg-C bonds. One hole distributes 
entirely over the 2p orbital associated to the long bond 
but the second spreads evenly over the remaining orbitals 
again giving a half- metallic ground state (FigJTJ)). ASIC 
changes the Jahn- Teller distortion into two short (2.15 A) 
and four long (2.18 A) bonds, basically reversing the 
LSDA geometry. In addition ASIC correctly describes 
the level occupation giving a filled singlet (shorter bond) 
separated by about 2.5 eV from an empty doublet. 

The triple acceptor Bo presents a somewhat less crit- 
ical situation. Now the B 2p impurity band is half- 
filled and positioned completely within the MgO band- 
gap. The strong Hund's interaction produces a magnetic 
ground state with a magnetic moment of 3 /ie/B for both 
LSDA and ASIC. For such an occupation the Jahn- Teller 
mechanism does not operate and the LSDA and ASIC 
bandstructures are qualitatively similar. 

When we look at the magnetic coupling we find for 
both N and C a situation very similar to that of Al in 
Si02, i.e. a large magnetic coupling in LSDA and prac- 
tically no coupling in ASIC. The situation is rather dif- 
ferent for B where even LSDA does not predict magnetic 
coupling. Does this mean that ferromagnetism cannot 
exist in these compounds? The answer to this question 
is not simple and essentially relates to the behaviour of 
the impurities with additional doping. For instance by 
adding a fractional electronic charge to Co one should 
find a situation intermediate between that on Co and 
Nq. Such a charge however feels a weaker nuclear poten- 



tial and certainly localizes less. Although it is likely that 
it will occupy one of the empty Co P orbitals, leading 
once again to a half- metallic band-structure, the quan- 
titative details of its exact localization are not easy to 
predict. Since the ferromagnetism is the results of the 
subtle interplay between localization, necessary for the 
magnetic moment formation, and de-localization, neces- 
sary for the magnetic coupling, we believe that a more 
accurate electronic structure method capable of taking 
all these factors into account is needed. 

We finally move our attention to Ga vacancy (Vca) 
and of Zn substitutional for Ga (ZnQ a ) in GaN. The 
first is again a triple acceptor and it was recently pro- 
posed as possible source of d° ferromagnetism 1 —, while 
the second introduces only one hole, which is predicted 
non-magnetic by LSDA. Our results, obtained with a 64 
atoms supercell are summarized again in table U where 
we can observe that LSDA and ASIC agree on the mag- 
netic ground state of Vca, but give respectively a non- 
magnetic and a magnetic 5=1/2 state for ZnQa- The 
non-magnetic state found in LSDA is resistant to super- 
cell size and persists to cells as large as 256 atoms, which 
are already large enough to decouple the impurities lo- 
cated in the periodic mirror cells. 

Let us consider first the case of Vca- Here the situation 
is somehow similar to that of MgO:N, where the magnetic 
state is reproduced already at the level of LSDA, but 
the spatial distribution of the moment changes radically 
when SIC is considered. In LSDA, in fact, the four N 
surrounding the vacancy move outwards so that the dis- 
tance with the vacant site is 2.1 A, about 8% longer than 
that between Ga and N. Such displacement is isotropic 
and the three holes spread uniformly over the four avail- 
able bond directions of the wurtzite lattice (see figure [3]). 
Thus the material turns out to be half metallio^. The 
ASIC picture is however rather different. Also in this 
case the SIC allows the hole to be localized. The ASIC 
relaxation leads to two bonds of 2. 14 A and two longer 
of bonds of slightly different lengths. The longest of the 
four is fully filled and the three holes localize on the re- 
maining three (Fig. [3]). Such a situation reminds that of 
Alsi and it is a further demonstration of the effects of 
self-trapping. Also for GaNiVca self-trapping massively 
suppresses the magnetic coupling energy which reduces 
from 145 meV (antiferromagnetic) in LSDA (for two Vc a 
placed at 6. 3 A from each other), to 1 meV for ASIC. 

The situation for the single acceptor Znc a is dramati- 
cally different with LSDA giving a non-magnetic ground 
state with the hole evenly distributed over the four al- 
most tetrahedrally coordinated bonds. These are just 
slightly expanded to about 2 A. ASIC again captures 
the Jahn- Teller distortion: three Zn-N bonds re- adjust 
to a distance of 1.96 A, while the forth expands consid- 
erably to 2.4 A. The hole now completely localizes over 
the longest bond, and due to the enhanced degree of lo- 
calization Hund's coupling spin splits the energy level 
in forming a spin 1/2 ground state. The ASIC electronic 
structure then is that of Znc a being a paramagnetic deep 





FIG. 3: (Color on line) Isosurface of the 3 holes (the magne- 
tization) associated to VGa in GaN as calculated with LSDA 
(left) and ASIC (right). In LSDA the holes are distributed 
uniformly over all the four N ions surrounding the vacancy, 
while they localize around the longer three bonds in ASIC. 




trap instead of a shallow non-magnetic acceptor. 



In conclusion we have demonstrated that the problem 
of predicting the magnetic ground state of p-type impu- 
rities is an extremely delicate one. In general standard 
LSDA is unable to capture Jahn- Teller distortion and sys- 
tematically underestimates the electron localization. As 
a consequence Hund's coupling is also underestimated. In 
the best case scenario a magnetic ground state is still pre- 
dicted, but the spatial distribution of the magnetic mo- 
ment is inaccurate. This typically leads to a metallic (via 
impurity band) ground state, strong exchange coupling 
between the magnetic impurities, and the expectation of 
room temperature ferromagnetism. The removal of the 
self-interaction error results in Jahn- Teller distortion and 
the creation of insulating ground state with little to no 
magnetic coupling between the impurities. A more dras- 
tic situation is however encountered when the erroneous 
LSDA derealization strongly suppresses Hund's coupling 
and produces a non-magnetic ground state. Also in this 
case ASIC is able to fix the problem and predicts non- 
vanishing magnetic moment localized at distorted bonds. 



FIG. 4: (Color on line) Isosurface of the hole (the magneti- 
zation) associated to Znc a in GaN as calculated with LSDA 
(left) and ASIC (right). In LSDA the hole is unpolarized and 
distributed uniformly over all the four N ions surrounding the 
vacancy, while it spin-splits and localizes around the longer 
bond in ASIC. 
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